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We investigate the behavior of a one-dimensional diatomic fluid under a shock wave excitation. 
We find that the properties of the resulting shock wave are in striking contrast with those predicted 
by hydrodynamic and kinetic approaches, e.g., the hydrodynamic profiles relax algebraically toward 
their equilibrium values. Deviations from local thermodynamic equilibrium are persistent, decaying 
as a power law of the distance to the shock layer. Non-equipartition is observed infinitely far from 
the shock wave, and the velocity-distribution moments exhibit multiscaling. These results question 
the validity of simple hydrodynamic theories to understand collective behavior in Id fluids. 

PACS numbers: 



Most one-dimensional (Id) systems -both quantum 
and classical- display anomalous collective properties. 
Generally speaking, the dimensional constrain inherent 
in these systems affects strongly the propagation of fluc- 
tuations. These fluctuations, in turn, control the sys- 
tem cooperative behavior, so their abnormal propagation 
leads to the observed anomalies. This observation under- 
lies most surprising results found in Id [l|,|2j]. 

Particularly interestingis the nonequilibrium behavior 
of Id classical fluids [iHSSEHBSHlIi An 
essential feature of these systems is that the particle se- 
quence remains unchanged during their time evolution. 
This important trait, also known as single-file constrain, 
is directly responsible for the anomalous propagation of 
fluctuations. For instance, in order for a particle in a Id 
fluid to move appreciably, a coherent, correlated motion 
of many particles is needed. This strongly suppresses dif- 
fusion: in Id (stochastic) fluids the mean square displace- 
ment grows now as (Ax 2 )sf = 2Ft 1 / 2 , much slower than 
usual Fickian diffusion (~ t) |3J • This single-file diffusion, 
which turns out to be relevant in many fields ranging 
from membrane biophysics to carbon nanotubes, zeolites, 
DNA, polymers and nanofluids, etc., has been recently 
confirmed in experiments with colloids 0]. Common to 
all these systems is the presence of confining structures 
(e.g. narrow channels) inducing quasi-ld (single-file) mo- 
tion. 

In the same way anomalous diffusion reflects the strong 
spatial correlations in Id fluids, energy transport in Id re- 
veals the presence of long-range velocity and current cor- 
relations. The simplest nonequilibrium situation where 
this is put forward is Fourier's law. When the Id fluid is 
subject to a small temperature gradient VT, an energy 
flux J = — kVT appears from the hot to the cold reser- 
voir. It has been found that the heat conductivity k has 
unusual large values for Id fluids. In particular, it is ac- 
tually believed (though some controversy remains 0,0) 
that any classical Id fluid with momentum-conserving 
interactions and non-zero total pressure should exhibit a 



divergent k in the thermodynamic limit, i.e. k oc L a as 
the system size L goes to oo. Although there exist ana- 
lytical and numerical evidences supporting this result in 
both Id fluids and crystals, there is still no agreement 
on the theoretical framework and the exponent a > 0: 
fluctuating hydrodynamics predicts a — 1/3 6], while 
mode-coupling theories and Boltzmann equation result 
in a = 2/5J3, , and numerical results are still incon- 
clusive 0, 0] (see however [j| for new strong evidences 
supporting a — 1/3). In this way, the behavior of Id 
fluids far from equilibrium and the theoretical approach 
suitable for their description are still open problems. 

In this Letter we study the response of a Id model fluid 
to a nonequilibrium excitation. In particular, we study 
the piston problem, a textbook example where we learn 
the notion of shock wave • Our results show that the 
anomalous nonequilibrium response of Id fluids cannot 
be described within a classical hydrodynamic approach. 
In particular, the two basic hypothesis underlying any 
hydrodynamic theory are violated in our model, namely: 
(i) perturbations away from local thermodynamic equi- 
librium decay algebraically, persisting for a long-time, 
and (ii) there are additional slowly-changing observables 
apart from the density, velocity, and energy fields. 

Our model is the diatomic gas of hard-point particles 
CH EJ - In a line of length L we introduce N = L hard- 
core point particles interacting via elastic collisions. The 
particles have alternating masses, m,2i—i = m = 1 and 
Tn>2i — M, i 6 [1, j]. In addition, we include a moving 
piston of infinite mass, which starts moving from x — 
at t = with constant velocity V = 1. All other parti- 
cles are initially at rest (cold gas)[ll|, and start with ran- 
dom positions, provided that the alternating order holds. 
All simulations here reported correspond to systems with 
N = 2 x 10 4 particles (other sizes were also simulated to 
check our results against finite-size effects), and we aver- 
age our results over 10 4 — 5 x 10 4 different realizations of 
the initial state. 

The propagation of the piston into the cold fluid ere- 
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FIG. 1: (color online) Number density profiles for different 
mass ratios. From bottom to top, n = 0.9, 0.8, 0.7, 0.618, 0.5. 
For fi < 0.5 profiles have been shifted 1 unit in the vertical 
axis. From top to bottom, fi = 0.4, 0.3, 0.2, 0.1, 0.01. 



ates a shock wave traveling with velocity W > V. At 
any time t, the flow field consists of an aftershock re- 
gion Vt < x < Wt of excited gas attached to the pis- 
ton, and a pre-shock region x > Wt of undisturbed gas. 
The position of the shock wave can be identified with 
the leading particle, trivially defined at any t as the last 
particle disturbed from its initial, motionless state. For 
equal masses, m = M, the shock structure is particularly 
simple. In this case the system is equivalent to a Id gas 
of non-interacting particles: the leading particle moves 
with velocity 2V, and the aftershock fluid is composed 
by particles with velocity and 2V, so the average is 
V. However, this limit is pathological, since the system 
is now non-ergodic, and lacks any relaxation mechanism 
in velocity space capable of driving the fluid toward (lo- 
cal) equilibrium. We may restore relaxation by letting 
M > to. In this case one expects to recover far behind 
the shock wave an equilibrium state with number density 
n, average flow velocity V and temperature T. Conser- 
vation of particles, total momentum and total energy, or 
equivalently the constancy of the corresponding fluxes 
(Rankine-Hugoniot conditions), thus imply. [f2T| 

n = 2n , W = 2V, T=-(m + M)V 2 , (1) 

where n is the density of the undisturbed gas. 

The randomness of the initial positions implies that 
the location and velocity of the shock wave (equivalently, 
the leader) are stochastic variables when consider within 
the ensemble of initial configurations. While the aver- 
age shock position grows as (xi) = Wt, the variance 
of the shock location is measured to grow diffusively, 
(Ax|) = 2Dt, with D an effective mobility (see |l3j 
for similar shock fluctuations). To prevent these fluc- 
tuations from blurring the shock wave structure^), we 
measure the flow profiles, characterizing the transition 
from the undisturbed gas toward the asymptotic equi- 
librium state in a reference frame moving with the 
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FIG. 2: (color online) Log-log scaling plot of An'(x) and 
Av'(x). Many different /j, £ [0.01,0.9] are included on each 
plot. The scaling factors have been selected in each case 
to fix the amplitude of the power-law tail to 1. 



leading particle. Moreover, profiles quickly converge to a 
steady shape after a short transient, thus allowing time 
averaging. Fig. ^ shows the normalized number density 
profiles measured in this way. Normalized observables 
are defined as 4>'{x) = (4>(x) — 0_)/(0+ — (/>_), where 

4>± = 4>{ x — * ±oo). 

Most interesting in Fig. ^ is the dependence of pro- 
files on mass ratio, fi — m/M. For fi close to 1, the 
profile changes smoothly; in particular, the shock width 
(the length scale on which profiles change appreciably) is 
much larger than the mean free path. However, as \x de- 
creases, the shock becomes steeper. This tendency lasts 
up to fi w 0.5. At this point a density overshoot devel- 
ops. Further decreasing /i results in a smoothening of the 
profile, with the exception of the initial overshoot, which 
sharpens as fi decreases. The shock width exhibits then 
a characteristic non-monotonous /i-dependence, with a 
minimum around fi ps 0.5. Notice that [i ~ 0.5 was pre- 
viously reported to yield the fastest relaxation times for 
this model .5]- Similar behavior is observed in velocity 
and energy profiles (not shown). Interestingly enough, 
both density and velocity overshoots are not observed in 
hard-sphere fluids for d > 1. 
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FIG. 3: (color online) P e (v) and Pf'iv), k = m, M for differ- 
ent mass ratios. 
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FIG. 4: (color online) Log-log scaling plot of 7m (a;) (main) 
and |£ m (a;)| (inset) versus the distance to the shock. Defini- 
tion of amplitudes and values of fi as in Fig. 



For large x, where x is the distance from the shock 
layer, one expects the fluid to relax toward the equilib- 
rium state (Q. To characterize this relaxation, we now 
measure the excess flow fields as a function of x, de- 
fined as A<p'(x) = 1 - (j)'{x), with <j>' = ri, v', E' . This 
is shown in Fig. [21 It is remarkable that these observ- 
ables relax toward the asymptotic state JD as a power 
law, Acp'(x) ~ x - ' 3 *. In particular, we find /3 n ~ f , and 
/3 V = Pe ~ 0.95 7^ /3 n (energy curves are not shown). 
These algebraic tails do not depend on /u, as demon- 
strated by the scaling in Fig. although the observed 
transient behavior is large for [i close to and 1, and min- 
imal for fi rj 0.5. The power laws are a reflection of the 
large correlations present in Id fluids, and are in striking 
contrast with the relaxation predicted by hydrodynam- 
ics and kinetic approaches, which is always exponential 
[ill Il5l]. In fact, this slow algebraic relaxation is related 
to the presence of long- wavelength hydrodynamic fluctu- 
ations |l8j . which control the fluid's transport properties. 
We discuss further this connection below. 

We now focus on the shock wave propagation. Fig. 
shows the shock stationary velocity distribution, Pi{v) = 
P{x — 0,v), for several /i's. The average propagation 
velocity is W = 2 in all cases, as expected. However, 
a multi-peaked structure in Pe(v) emerges for fi < 0.5, 
which coincides with the appearance of profile over- 
shoots, see Fig. This exemplifies how far the system 
is from local thermodynamic equilibrium near the shock 



layer. We may split Pe(v) = Pj (v) + P^ vl> (v), where 
Pt\v) is the probability density function (pdf) for find- 
ing a particle of mass k = m, M and velocity v leading 
the shock (or equivalently, the pdf P( k )(a; — 0, v) for hav- 
ing a particle of mass k — m, M at x = with velocity v 
in a reference frame moving with the shock wave). These 
distributions, plotted also in Fig. show that the shock 
wave propagates faster (slower) whenever a light (heavy) 
particle is leading. Also remarkable is the incipient sin- 
gularity in Pi(v) for small v as fi —> 0. This, together 



>(M), 



with the singular overshoot in the profiles for small fi 
(Fig. nj, makes the limits \i — ► and fi — * 1 essentially 
different, contrary to recent claims 

Next, we turn our attention to the emergence of local 
thermodynamic equilibrium (LTE) behind the shock. For 
that, we measure the skewness 7k(a0 and kurtosis excess 
|£k(x)| of the local velocity distribution for each species, 
P( k \x,v), k — m,M. They are defined as 7k(#) = 



[4 k) (4 k) y 3/2 ]i x ) and 



}(x), 

where v^(x] is the i— th central moment of P^ (x,v). 
As the system reaches LTE behind the shock, both 7k (x) 
and |£k(aO| must go to zero. Fig. 0] shows the actual 
relaxation. Surprisingly we find that both quantities de- 
cay also as slow power laws, (j>k(x) ~ a; - ' 3 *, where now 
(3 1 w 0.95 and /3t ~ 0.4 (exponents are independent of 
the species k = m,M). This proves that perturbations 
away from LTE are persistent in Id, and that higher- 
order fields other than the hydrodynamic ones (namely, 
7k (x) and |^k(^)|) change slowly and therefore are rele- 
vant. Also interesting is the fact that /3 1 > This ob- 
servation, combined with the other exponents f3 n ,v,E, im- 
plies that the distributions P^ (x, v) exhibit multiscal- 
ing. This emerging picture is intriguing and exemplifies 
the anomalous behavior of Id fluids. Another interesting 
observable is the total energy stored in light and heavy 
particles at a given time t, e^(t), k — m,M. Energy 
equipartition in the asymptotic equilibrium state implies 
that ek(i) — * ;j(3fc + k)V 2 = ek(oo) as t — > 00, where 
k is the mass other than k. The inset to Fig. [5] shows 
the excess energies Aet(i) = £k(oo) — ek(t), k = m,M, 
for a particular /i. They also decay as a power law, 
as expected from the algebraic tail of AE'(x). How- 
ever, what is remarkable is the different amplitudes: light 
particles are always closer to the asymptotic equiparti- 
tion energy than heavy particles (ie. lig ht particles are 
more energetic than heavy ones 0, Il5l fl6l fl7j|). Wc 
may characterize this non-equipartition of energy with 
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FIG. 5: (color online) Main: amplitude ratio X(/j.) (see main 
text for definition) as a function of [i. Inset: Log-log plot of 
Aek(t), k — m,M, vs. time for the particular case fi = 0.8. 
The dashed line has slope 0.95 « (3e- 



suggests that the anomalous collective properties of Id 
fluids cannot be described within a simple hydrodynamic 
approach. This is particularly important in the study of 
transport phenomena in these systems, where hydrody- 
namic equations have been used to predict an anomalous 
thermal conductivity, for instance |fg. A sound theoret- 
ical approach to collective phenomena in Id fluids must 
therefore go beyond hydrodynamics. A good candidate is 
generalized kinetic theory. However, taking into account 
the strong inter-particle correlations present in Id fluids 
in such a microscopic theory is a challenge that deserves 
further attention. 

The author thanks P.L. Krapivsky for the original idea 
of this study, and T. Antal, P.L. Garrido, P.L. Krapivsky, 
and S. Redner for useful discussions. Financial support 
from the Spanish MEC is also acknowledged. 



the ratio A(/z) = Ae m (t) / AemW > which converges to a 
/i-dependent constant for long enough times. This is 
shown in Fig. \5\ where it is evident that energy non- 
equipartition becomes more important as /i decreases. 

A full theoretical explanation of the observed behavior 
seems challenging. Both hydrodynamic [10( and (Boltz- 
mann) kinetic 15] approaches do not account for the 
measured shock structure and the slow, power-law relax- 
ation of hydrodynamic observables. These algebraic tails 
result from long-wavelength hydrodynamic fluctuations, 
a main source of non-locality and long-range correlations 
in fluids. In fact, for 3d fluids reflection and refraction 
of hydrodynamic fluctuations on a shock wave give rise 
to a slow, universal power-law relaxation of all hydrody- 
namic fields, Acf>(x) ~ x~ 3 / 2 , cj> — n,v,E 18]. Although 
somewhat similar in spirit, Id fluids exhibit a much more 
intricate behavior. The algebraic tails for the hydrody- 
namic fields are not universal, with exponents /3 n rs |, 
and /? v = Pe ~ 0.95 ^ (3 n . In addition, we observe non- 
equipartition of energy and multiscaling behavior of one- 
particle distribution functions. Most importantly, devi- 
ations away from local thermodynamic equilibrium are 
persistent, as characterized by the slow power-law de- 
cay of the skewness and kurtosis fields, with exponents 
/3 7 k 0.95 and f3^ « 0.4, respectively. 

Any hydrodynamic description of collective behavior 
relies on two basic hypothesis jg: (i) the system must 
reach local thermodynamic equilibrium rapidly (i.e. ex- 
ponentially fast), and (ii) the only slowly-evolving ob- 
servables are the density, velocity and energy fields. The 
slow algebraic relaxation of Jk(x) and |£k(x)| proves that 
both hypothesis are violated in our Id model fluid, which 
on the other hand contains the main ingredients charac- 
terizing a broad class of Id systems (i.e. dimensional con- 
strain, and energy and momentum conservation). This 
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